addpath('../lib');
addpath('../mockData');

disp('initial computations');


niter = 130;

t = 0:0.3:10; %laikas
lam = 100:10:700; %bangos ilgiai

disp(size(t));
disp(size(lam));

spec = MockWGaussian(lam);
conc = MockHA2B2C(t);


duomenys=spec*conc;
R = 3;
%% standart NMF
disp('Computing NMF');

[W1 H1] = TNMF(duomenys, R, niter);
disp('Finished');
